Infrared observations of a mid-L dwarf with strong Ha emission 
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ABSTRACT 

We present Spitzer/IRAC observations of the L5 dwarf, 2MASSI J1315309- 
264951 (2M1315). This ultracool dwarf is known to display strong emission in 
the Ha line. The SED for this object does not show any IR excess, that would 
indicate the presence of an accretion disk. Although the IRAC colors for 2M1315 
are consistent with other L dwarfs, they seem to be redder by ~0.1 mag compared 
to the other L5 dwarfs, and more like the late-type L dwarfs. The existing six 
epochs of spectroscopy suggest that the emission in Ha is not persistent, but 
shows long-term variability between a flare value of ~100 A and a quiescent 
value of ~25 A. Chromospheric activity seems to be the most likely cause, which 
is also indicated by the detection of Na I D lines in emission (Fuhrmeister et 
al.). We have measured a proper motion of 0.79"± 0.06"/yr, that corresponds 
to a tangential velocity of ~81 km/s, at a distance of ~22 pc. The high Vt an for 
this object suggests an old age. Evolutionary models indicate lower limits of 3.3 
Gyr (Chabrier et al.) or 1.4 Gyr (Burrows et al.) for this magnetically active L 
dwarf to be a stable hydrogen-burning star. There seem to be no observational 
differences between this old L dwarf that has Ha emission and the other field L 
dwarfs without it. 

Subject headings: stars: brown dwarfs - stars: activity - stars: chromospheres 



1. Introduction 

Emission in Ha is commonly seen among M dwarfs, with the activity levels increasing 
towards later types due to longer spin-down timescales, reaching a peak at M7. Beyond 
spectral type M7, there is a drop off in Ha activity, with a rapid decline seen for L dwarfs 
(Gizis et al. 2000). This suggests a decrease in the generation and dissipation of magnetic 
fields for late-M and L dwarfs. One reason could be the increasing atmospheric neutrality in 
ultracool dwarfs that decouples the field from atmospheric motions, and can hamper their 
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dissipation, even if strong fields are generated (Mohanty et al. 2002). This decoupling can 
cause suppression of chromospheric and coronal heating, resulting in a decrease in emission 
in the Ha line and X-rays, despite the presence of magnetic fields. Indeed, magnetic field 
strengths of 0.1-lkG have been measured for some late-M and L dwarfs (e.g. Berger 2006) 
indicating that the drop in Ha emission seen for ultracool dwarfs must be due to atmospheric 
neutrality rather than a fall in magnetic flux production (Reiners & Basri 2006). 

On the contrary, there exists a small population of late-M and L dwarfs, such as 
PC0025+0447 (M9.5; Schneider et al. 1991), 2MASSW J0149090+295613 (M9.5; Liebert 
et al. 1999) and 2MASS J01443536-0716142 (L5; Liebert et al. 2003), that exhibit strong 
flaring in the Ha emission line. This flaring is occasional and occurs only on the order of 
~1% of the time (Liebert et al. 2003). While different explanations have been suggested 
for these rare displays of strong Ha emission, the mechanism that empowers such flares is 
still unknown. We report here on Spitzer (Werner et al. 2004) IRAC (Fazio et al. 2004) ob- 
servations of the mid-type L dwarf 2MASSI J1315309-264951 (hereafter, 2M1315). 2M1315 
was found to show strong emission in Ha by Hall (2002a), that decreased by a factor of 2 
half a year later (from an equivalent width (EW) of 121±3lA to 25±10A). Hall classified 
this object as an L3 dwarf, and concluded that the variable Ha emission is either due to 
chromospheric activity, or by accretion in a binary system. Gizis (2002) independently found 
2M1315 and classified it as an L5 dwarf. However, he found this object to be a persistent 
strong Ha emitter (EW =97±10A) from 6 consecutive 600s exposures on one night, and 2 
consecutive 900s exposures the following night. On the basis of an absence of low surface 
gravity signatures in its spectrum, Gizis concluded that this object is not a member of TWA, 
and may be a relatively old L dwarf. This was confirmed by Hall (2002b) who found 2M1315 
to have a high proper motion (0.71"± 0.07"/yr), that corresponds to a V tan of ~76 km/s at 
d~23 pc. Hall (2002b) again found a large EW of 124 A, and thus concluded that 2M1315 is 
a persistent strong Ha emitter. Later, Fuhrmeister et al. (2005) reported the weak detection 
of Na I D lines in emission, and an Ha EW of 24.1 A, suggesting that 2M1315 was in a qui- 
escent state since Hall's measurement. Recently, Barrado y Navascues (2006) has reported 
an average EW of 153±26 A from seven consecutive spectra obtained in March 2003. 

The likely explanations for strong Ha emission can be narrowed down to (i) presence 
of an unusually strong magnetic field, (ii) accretion in a closely interacting brown dwarf 
binary system, (iii) presence of an accretion disk. We have obtained IRAC observations of 
2M1315 in order to check for IR excess, that would indicate the presence of a disk around 
it. In Section [21 we present our observations. Section [3] gives a discussion of the previous six 
epochs of spectroscopy, and the likely causes of strong Ha emission seen in this L dwarf. 
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2. Observations 

Observations in the IRAC bands were obtained in February, 2006, using the full array 
mode. A frame time of 12 sec was requested, along with a dither pattern of five positions 
Gaussian, which allows the removal of instrumental artifacts. Aperture photometry was 
performed using the task PHOT under the IRAF APPHOT package. An aperture radius of 
5 pixels and a background annulus of 10-20 pixels was used. The zero point fluxes of 280.9, 
179.7, 115 and 64.1 Jy and aperture corrections of 1.05, 1.05, 1.058 and 1.068 were used 
for IRAC channels 1 through 4, respectively. The errors in magnitudes are between ~0.01 
and 0.02 mag, with an additional IRAC calibration uncertainty of ~2%. Table 1 lists the 
observed fluxes in the IRAC bands. 



3. Discussion 

In order to determine the extent of IR excess, we have looked at different models for very 
low-mass stars so as to obtain a good fit to the photospheric emission from 2M1315. The 
models by Allard et al. (2001) are two limiting cases of dust treatment in the atmosphere. 
The DUSTY models include not only grain formation, but also grain opacities. As first 
shown by Tsuji et al. (1996), dust opacities significantly affect the emergent flux from an 
ultracool dwarf, as they deplete the gas phase of a number of molecular species, such as, TiO, 
VO, FeH, CaH, MgH, and of refractory elements, such as Al, Ca, Ti and Fe. The resulting 
spectral distribution thus appears redder with weaker molecular opacities, as compared to 
the models that do not consider grain opacities. The COND models by Allard et al. (2001), 
do consider grain formation but do not take into account the grain opacities, resulting in a 
spectral distribution that is transparent blueward of l.O^m, and very similar to the dust- 
free case (the NextGen models; Hauschildt et al. 1999a,b). As discussed in Chabrier et al. 
(2000; hereafter COO), in the DUSTY models, turbulent mixing is very efficient in bringing 
material upward to the region of condensation and maintaining small grain layers, which 
would otherwise gravitationally settle down below the photosphere. In the COND model, 
however, the turbulent mixing cannot balance gravitational settling, which results in the 
sedimentation of all of the condensed species below the photosphere. Intermediate between 
the two extremes of DUSTY and COND are the SETTLE models by Allard et al. (2003). 
These models assume a sedimentation efficiency higher than the DUSTY but lower than the 
COND models, i.e. when grains have partially settled below the photosphere. 

A spectral type of L5 (Gizis 2002) suggests an effective temperature of ~1700 K (Kirk- 
patrick 2005) for 2M1315. In Fig. la, we have used the DUSTY, SETTLE and COND 
models to get a good fit to the photospheric emission from 2M1315. The COND models 
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result in bluer colors for decreasing T e ff, and more flux shortward of 1/zm, as compared to 
the DUSTY ones, while the SETTLE models show fluxes intermediate between the two. 
Fig. la shows the observed near-IR flux to be over-estimated by the COND and SETTLE 
models for 2M1315, whereas the DUSTY one provides a good fit. Allard et al. (2001) have 
discussed the affects of gravity on the DUSTY models (see their Figs. 16 and 17). Due to the 
higher pressures and densities of the high-gravity models, the optical to red spectral region 
is significantly veiled by dust opacities. The near-IR water vapor bands are also affected. A 
log g=5.0 DUSTY model underestimates the observed flux in the J-band for 2M1315, as can 
be seen from Fig. lb, whereas log g=4.0 provides a good fit. We note that a lower value of 
log g does not imply that 2M1315 is a relatively young star. The temperature 1700 K lies at 
the cut-off between the DUSTY and the COND models, suggesting that the models may be 
inconsistent at this temperature. The observed (J-K) colors for the reddest L dwarfs show 
a saturation at a value of ~1.9, corresponding to a T e ff ^1700-1800 K (see Fig. 6 in COO). 
2M1315 lies in this region, with a (J-K) color of 1.7 and an absolute K magnitude of 11.8 
(Fig. lc). Due to the continued influence of dust opacities, the DUSTY models continue to 
redden to colors (J-K)>2, and thus fail to reproduce the colors for mid- to late-type L dwarfs 
(M K >11.5, J-K >1.5). COO have explained this to be caused by the gravitational settling 
of some grains below the photosphere. Fig. lc shows that while the SETTLE models can 
reproduce the turn to the blue of (J-K) colors that is seen from mid-L to T-types in color- 
magnitude diagrams, they still fall short of reproducing the colors of the reddest L dwarfs. 
The cloudy models by Marley et al. (2002) which, similar to the SETTLE case, assume a 
sedimentation efficiency intermediate between the DUSTY and COND models, have been 
able to reproduce the (J-K) colors for the latest type L dwarfs. 

The SED for 2M1315 (Fig. la) does not show any significant excess emission in any 
of the IRAC bands. The 3.6 to 8 /im slope for 2M1315 is ~-2.5. According to the criteria 
defined in Lada et al. (2006), this makes it a disk less object, consistent with our finding. 
The main spectral features in the IRAC bands are due to molecular absorption of CH 4 near 
3.3 and 7.8 fim, and H 2 at 5.5 and 6.5/im (Noll et al. 2000, Roelling et al. 2004, Cushing 
et al. 2006). The absorption in these bands strengthens towards later types. The dip seen 
in 2M1315's SED at 4.5/im could be due to the CO band at 4.7/im, as seen in the T dwarf 
GL 229B (Oppenheimer et al. 1998). 

In Fig. 2, we have compared the IRAC colors of 2M1315 with other M, L and T dwarfs 
from Patten et al. (2006). Although the colors for 2M1315 are consistent with other L 
dwarfs, they seem to be redder by ~0.1 mag compared to the other L5 dwarfs, and more 
like the late-type L dwarfs. More CO absorption near 4.5 /im relative to H 2 near 5.5 /im 
can explain the redder [4.5]-[5.8] color. A redder [3.6]- [8] color suggests more absorption of 
methane near 3.3//m, compared to 7.8/mi. Could a (slightly) redder color suggest a higher 
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surface gravity for 2M1315? To answer this, we have checked for any variations in Spitzer 
colors with surface gravity or tangential velocity. In Fig. 2, L dwarfs from Patten et al. for 
which Li EW has been measured are denoted by filled diamonds. For the others, Li EW 
could not be found in the literature. We have also determined the tangential velocities for 
M and L dwarfs in Patten's sample by searching for their proper motion and parallaxes in 
the SIMBAD database. We could only find proper motions for 24 objects; out of these, six 
have V tan between 50 and 80 km/s, and are denoted by filled circles. Fig. 2 shows that 
the IRAC colors for M and L dwarfs are not very sensitive to surface gravity, or tangential 
velocities. Though the number of data points used for comparison is very small, the lack of 
any concentration of the higher velocity (older) objects near 2M1315, or that of the objects 
with lithium (lower surface gravity, younger) near bluer IRAC colors suggests that these 
colors are not very useful as age indicators. 

From the 2MASS and IRAC images, we have determined the proper motion for 2M1315 
to be 0.79"± 0.06"/yr (/i a = -0.73"/yr, fi s = -0.31"/yr). Earlier, Hall (2002b) had found a 
proper motion of 0.71"±0.07"/yr, assuming that the nearby bright star, USNO J131531.230- 
264953.01, that lies at a separation of 7.14" ± 0.05" (epoch 2002.685) has negligible proper 
motion. Comparing the 2MASS and IRAC images shows that the position of the USNO star 
has not changed, and thus confirms that it is 2M1315 that exhibits the high proper motion. 
A distance of 21.7 pc (obtained from the relation in Dahn et al. 2002) results in a tangential 
velocity of 81.3 km/s. The proper motions for ten known TWA members range between 0.06 
and 0.13"/yr, resulting in V tan between 13 and 32 km/s (using distances from Reid 2003). 
This rules out a TWA membership for 2M1315, and suggests an old age. That 2M1315 is 
relatively old and massive was suggested earlier by Gizis (2002) and Hall (2002b), based 
on the absence of low surface gravity signatures and its high proper motion, respectively. 
Fig. 3 shows the age isochrones for masses O.O9M (top) to O.O2M (bottom), from COO 
DUSTY (bold lines), COND (long-dashed lines), and Burrows et al. (1997; hereafter B97) 
(short-dashed lines) models. The B97 models are updated versions of their 1995 models, and 
show a change of less than 10% in luminosity at a given age and mass. COO have discussed 
that the number of grain species in the B97 model is very small, and the grain opacities are 
treated as frequency- independent, resulting in a gray atmosphere structure above 1300K. COO 
models consist of a significant number of grain species, the opacities of which are frequency- 
dependent, resulting in nongray atmosphere for 900 < T e ff < 2800K. Fig. 3 suggests that 
the decrease in effective temperatures due to silicate formation between 10 8 and 10 9 yrs is 
steeper in B97 tracks, as compared to COO models. This could explain the lower effective 
temperatures seen for B97 tracks during the H-burning phase. At a given age and mass, 
the temperature difference between the two models can be as large as ~400K, and seems 
to be more pronounced for the higher masses. The H-burning limit for the DUSTY models 
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is 0.07M Q . COO had found that including dust opacities lowers the effective temperature 
and luminosity for a given mass. The difference in T e ff between the DUSTY and COND 
models is more obvious at lower temperatures, and could reach up to ~10%. The COND 
models thus have a higher H-burning limit at O.O72M . The Li-burning limit for a 1 Gyr 
object from DUSTY models is O.O55M . For higher masses, the Li abundance goes to zero. 
Hall (2002b) had suggested an age of > 2 Gyr for 2M1315. For an effective temperature of 
1700 K, 2M1315 would have to be younger than ~760 Myr (COO) or ~630 Myr (B97) to be 
a Li-burning brown dwarf. For it to be a H-burning star, the lower age limits are 3.3 Gyr 
(COO) or 1.4 Gyr (B97). The strong Ha emission shown by this likely H-burning star is in 
agreement with the results from Gizis (2000), who concluded that some active L dwarfs are 
drawn from an older, more massive population. He has also noted the presence of another 
population of L dwarfs like Kelu 1, that show both Ha emission and Li absorption. Given 
2M1315's old age, its emission can be related to having a higher surface gravity and a fusion 
energy source. 

There are now six reported measurements of Ha EW for 2M1315. These are listed in 
Table 2, along with the dates of observation and the spectral resolution. Hall (2002a) have 
discussed that their 2001 August EW measurement of 25 ± 10 A is an underestimate as 
that spectrum may have been contaminated by the flux from the nearby USNO star. Thus 
discarding this measurement, he concluded that 2M1315 shows persistent non-variable Ha 
emission of EW ~ 100 A, with log (L Ha /L bo i) ~ -4. However, the EW measurement of 24.1 
A by Fuhrmeister et al. (2005) is from a ~45,000 resolution spectrum, and indicates that 
the Ha emission seen for 2M1315 may not be persistent, but variable. The values in Table 2 
suggest a long-term variation in Ha emission, with an EW of ~25 and ~100 A defining the 
quiescent and the maximum states, respectively. We note that the Ha EW depends strongly 
on the resolution, as discussed in White & Basri (2003). At low resolution, the 6569 A TiO 
feature blends with the Ha line, resulting in an underestimate of the continuum redward of 
the Ha feature, and an overestimate of the Ha EW. The higher the resolution, the smaller 
the EW, with the changes being as large as ~30% in going from low to high resolution. 

Our Spitzer observations can now rule out the presence of an accretion disk around 
2M1315, which could have explained its strong Ha emission. An absence of IR excess was, 
however, expected, considering the conclusions made by Gizis (2002) and Hall (2002b) about 
its age. An age of a Gyr or more makes 2M1315 unlikely to be a T Tauri analog brown dwarf, 
and is old enough to sweep away all signs of a debris disk around it. 

The other possible explanation for Ha emission in 2M1315 could be the accreting brown 
dwarf binary scenario, discussed further in Burgasser et al. (2000). However, as Hall (2002a) 
suggested, this would require some variability in the accretion from the sustained Roche lobe 
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flow, that would explain the variability in Ha emission. 

The most likely explanation seems to be chromospheric activity. The long-term varia- 
tions in Ha suggest a slowly varying magnetic field, the strength of which though must be 
much stronger than that for other L and T dwarfs. This can be seen from Fig. 2 in Liebert 
et al. (2003), where even the quiescent emission of 2M1315 is higher by ~2 orders of mag- 
nitude than the upper limits for other L dwarfs at the same spectral subclass. Hall (2002a) 
had suggested that the variable Ha emission is more likely to be due to rapid variations, 
like flaring, than slow variations in the magnetic field strength, since this is rare in late-type 
dwarfs. However, flaring would require a shorter timescale for the variations to occur than 
that seen for 2M1315. Fuhrmeister et al. (2005) reported weak detection of Na I D lines 
in emission in 2M1315, the first such detection in an L dwarf. This again suggests that 
2M1315 is a chromospherically active object. Comparing with the other two very late-type 
field dwarfs, PC 0025+0447 and 2MASSI J1237392+652615, that show persistent strong Ha 
emission, Liebert et al. (2003) had concluded that these, along with 2M1315, could either 
be young (10-100 Myr) rare objects in the solar neighborhood, or a "special", yet unknown 
dynamo could be responsible for generating a stronger-than-average magnetic field in these 
objects. Such a possibility is discussed in Mullan & MacDonald (2001), who have considered 
several nonstandard additions to the "standard models" of main-sequence stars, in order to 
explain why the magnetically active and inactive stars form two distinct populations in a 
T e //-radius diagram. Their models suggest that the onset of a radiative core in an object of 
mass as low as O.15M requires T e ff to drop to a value of ~1500 K, and stronger magnetic 
fields with upper limits on strength of about 100 MG. In other words, a strong magnetic 
field could make a H-burning star look like a cool L dwarf. Mullan & MacDonald have also 
suggested that at a given T e ff, active stars have larger radii than inactive ones. 2M1315 
could be more luminous due to its stronger activity, if these nonstandard models are correct. 

We also consider the possibility that if the inclusion of dust opacities significantly influ- 
ences the emergent flux, it might also affect the chromospheric line emission. In their model 
chromospheres for mid- to late-type M dwarfs, Fuhrmeister et al. (2005) have compared the 
Ha line emission for the models with and without dust (see their Fig. 3). These authors 
have found that including dust opacities enhances most of the emission lines and lowers the 
continuum. The flux in the Ha line in their models increased by a factor of ~1.3 by including 
dust opacities. Hall (2002a) had found that the Ha flux in their 2001 March measurement 
(EW = 121 ± 31 A) for 2M1315 was a factor of 2 higher than the 2001 August measurement 
(EW = 25 ± 10 A). Thus dust opacities could affect the chromospheric line emission, and 
consequently the EW measurement. 

In summary, the Ha emission in 2M1315 is not persistent, but displays a long-term 
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variation between a flare value of ~100 A and a quiescent value of ~25 A. The absence 
of circumstellar material supports the chromospheric model. The magnetic field heating up 
the chromosphere must be unusually strong, but slowly varying. It is though not clear why 
this object has strong Ha emission from the chromosphere, while other old L dwarfs do not. 
2M1315 has a high proper motion (0.79"± 0.06"/yr) which results in a Vt an of ~81 km/s 
at a distance of ~22 pc. The high velocity suggests that an old age is likely, and rules out 
a TWA membership. At an effective temperature of 1700 K, evolutionary models indicate 
lower limits of 3.3 Gyr (COO) or 1.5 Gyr (B97) for this magnetically active brown dwarf to 
be a hydrogen-burning star. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by NASA through an award issued 
by JPL/Caltech. This work has made use of the SIMBAD database. 
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Fig. 1.— Top (a) DUSTY (blue), SETTLE (green) and COND (red) model fits for 2M1315. 
JHK magnitudes are from 2MASS, IRAC from this work. The errors in IRAC magnitudes 
are between 0.02 and 0.06 mag, and are smaller than the symbol size. Middle (b) DUSTY 
model fits for log g=4 (red) and log g=5 (green) for a T e ff of 1700 K. Bottom (c) Near-IR 
CMD comparing the three models: solid and long- dashed lines, DUSTY and COND models 
for 1 Gyr, respectively; short-dashed line, SETTLE models mapped onto the DUSTY 1 Gyr 
T e jj-radius relation (from Allard et al. 2003). Filled square, 2M1315; open squares, filled 
triangles and crosses, M, L and T dwarfs from Dahn et al. (2002), respectively 
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Fig. 2.— IRAC colors for 2M1315, along with M, L and T dwarfs from Patten et al. (2006). 
Filled square, 2M1315; open triangles, M dwarfs; open circles, L0-L5; open squares, L6-L8; 
plus signs, T dwarfs; filled diamonds, L dwarfs with Li EW measurements available; filled 
circles, M and L dwarfs with V tan between 50 and 80 km/s. 



-13- 



3000 



2800 



2600 



2400 



2200 



01 

I- 



2000 - 



1800 



1600 - 



1400 



1200 




-0.2 0.2 

log age (Gyr) 



0.4 0.6 



Fig. 3. — Evolutionary tracks for masses O.O9M (top) to O.O2M (bottom). Solid lines, 
DUSTY (COO); long-dashed lines, COND (COO); short-dashed lines, B97. 



-14- 



Table 1. IRAC Observations 



Star 


3.6 fim 


S/N 


4.5 fj.m 


S/N 


5.8 fim 


S/N 


8 fim 


S/N 




mjy 




mjy 




mjy 




mjy 




2M1315 


3.10±0.06 


60.86 


2.05±0.04 


60.19 


1.92±0.02 


119.53 


1.33±0.02 


57.32 



Note. — The 1-u errors do not include IRAC calibration uncertainty of ~2%. 



Table 2. Six epochs of spectroscopy for 2M1315 



Date of Observation 


Ha EW (A) 


Spectral Rcsolution(A) 


Ref 


2001 March 30 


121 ± 31 


8.6 


Hall (2002a) 


2001 May 1-2 


97±10 


2.8 


Gizis (2002) 


2001 August 15 


25 ± 10 


8.6 


Hall (2002a) 


2002 March 15 


24.1 


0.1 


Fuhrmeister et al. (2005) 


2002 September 9 


124 ± 55 


4.6 


Hall (2002b) 


2003 March 11 


153 ± 26 


2.7 


Barrado y Navascues (2006) 



